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The fact that the capacitance coefficients for a set of conductors are geometrical factors is de- 
rived in most electricity and magnetism textbooks. We present an alternative derivation based on 
Laplace's equation that is accessible for an intermediate course on electricity and magnetism. The 
properties of Laplace's equation permits to prove many properties of the capacitance matrix. Some 
examples are given to illustrate the usefulness of such properties. 



I. INTRODUCTION 

The fact that the capacitance is a geometrical factor 
is an important property in courses on electricity and 
magnetism. 1 ' 2 Derivations of this property are usually 
based on the principle of superposition^ and the Green 
function formalism^. Nevertheless, such derivations are 
not convenient for calculations. Alternative techniques 
to calculate the capacitance coefficients based on the 
Green function formalism^ and other methods^^ have 
been developed. 

In this paper we give a simple proof of the geometrical 
nature of the capacitance coefficients based on Laplace's 
equation. Our approach permits to demonstrate many 
properties of the capacitance matrix. The method is 
illustrated by reproducing some well known results, and 
applications in complex situations are suggested. 

II. CAPACITANCE COEFFICIENTS 

We consider a system of N internal conductors and an 
external conductor that encloses them. The potential on 
each internal conductor is denoted by tfi, i = 1, 2, . . . , N. 
The surface of the external conductor is denoted by 
Sjv+i, and its potential is denoted by tfN+i (see Fig.Q}. 
One reason to introduce the external conductor is that it 
provides a closed boundary to ensure the uniqueness of 
the solutions. In addition, many capacitors contain an 
enclosing conductor as for the case of spherical concen- 
tric shells. As we shall see, the case in which there is no 
external conductor can be obtained in the appropriate 
limit. 

The surface charge density a on an electrostatic con- 
ductor is given bji2 

o-i = e E • ih = -e V(f> ■ n* (i = 1, . . . ,N + 1), (1) 

where is an unit vector normal to the surface Si point- 
ing outward with respect to the conductor (see Fig. [T|) ; 
E and cj) denote the electrostatic field and potential re- 
spectively. The charge on each conductor is given by 

Qi = f o-i dS = -e f V0 ■ n., dS. (2) 
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FIG. 1: A system consisting of N internal conductors with 
conductor TV + 1 enclosing them. The normals m with i — 
1, . . . , jV+1 point outward with respect to the conductors and 
inward with respect to the volume Vs T (defined by the region 
in white). The surfaces Si with i = 1,...,7V are slightly 
bigger than the ones corresponding to the conductors. In 
contrast, the surface Sn+i is slightly smaller than the surface 
of the external conductor. 

The surface Si encloses the conductor i and is arbitrar- 
ily near and locally parallel to the real surface of the 
conductor (see Fig. [T])i2 We define the total surface St 
as 

St = Si + . . . + Sn + Sjv+i. (3) 

The volume Vs T defined by the surface St is the one 
delimited by the external surface Sn+i and the N inter- 
nal surfaces Si. The potential <fi in such a volume must 
satisfy Laplace's equation with the boundary conditions 

4>{S i ) = y i (i = l,...,N + l). (4) 

Because of the linearity of Laplace's equation, the solu- 
tion for <p can be parameterized as 

N+l 

< t>=Yl Vjfj, (5) 

where the fj are functions that satisfy Laplace's equation 
in the volume Vs T with the boundary conditions 

V 2 /, =0, = <%, (i,j = l,...,N + l). (6) 
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The solutions for fj ensure that is the solution of 
Laplace's equation with the boundary conditions in 
Eq. ((4]). The uniqueness theorem also ensures that the 
solution for each fj is unique (as is the solution for (f>). 
The boundary conditions ^ indicate that the fj func- 
tions depend only on the geometry. 

If we apply the gradient operator in Eq. (JSJ) and sub- 
stitute the result into Eq. ||2J), we obtain 



N+l 

<3i = E Gi Mi 



C l3 = -e <j> Vfj • n l dS, 



(7a) 
(7b) 



which shows that the Cij factors are exclusively geo- 
metric. The symmetry of the associated Cij matrix can 
be obtained by purely geometrical arguments. We start 
from the definition of Cy in Eq. (|7b|) and find 

Cij = -e / Vfj ■ m dS = s (s <f> fiWfj ■ i-ni) dS, (8) 

JSi JSt 

where we have used the fact that fa = 1 on the surface 
Si and zero on the other surfaces. From Gauss's theorem 
we obtain 

Cij=e [ V-(fiWfj)dV (9a) 

JVs T 

= e f [V/, • Vfj + f t V 2 f 3 } dV. (9b) 

JV.<!„ 



Because V 2 /, = in Vs T , it follows that 
Cij =e f V/, • V/, dV, 



(10) 



Equation (flTJ)) implies that is symmetric^ that is, 



Cij — Cji. 



(11) 



For certain configuration of conductors, consider two 
sets of charges and potentials {Qi,ipi} and </?■}. 
From Eqs. ([7]) and (jTTJ) we have that 

N + l N+l N+l 

E^i = E(E c ^ ( 12a ) 

i=l i=l j=l 

N+l N+l 

j=i i=i 

which implies that 

N+l N+l 

E Qifi = E QjVi' (1 3 ) 

i=l 3=1 

Equation p3p is known as the reciprocity theorem, 1 



When one or more of the N internal conductors has 
an empty cavity, is well known that there is no charge 
induced on the surface of the cavity^ 2 - (let us call it Si. c ). 
Consequently, although c is part of the surface of the 
conductor, such a surface can be excluded in the integra- 
tion in Eq. ([2]). In addition, we can check by uniqueness 
that fj = Sij in the volume of the cavity Vi t c so that 
V/j = in such a volume, and hence it can be excluded 
from the volume integral (|10[) . In conclusion neither S 1 ^ c 
nor Vi t c contribute in this case. 

The situation is different if there is another conductor 
in the cavity. In this case, the surface of the cavity con- 
tributes in Eq. ([5]). Similarly the volume between the 
cavity and the embedded conductor contributes in the 
volume integral (fTU)) . The arguments can be extended for 
successive embedding of conductors in cavities as shown 
by Fig. [5] or for conductors with several cavities. 




FIG. 2: Example of system in which there is a successive 
embedding of conductors. The volume Vs T corresponds to 
the region in white. The regions corresponding to empty 
cavities (and their associated surfaces and volumes) can be 
excluded without affecting the calculations. In this picture 
cavity A is empty and its surface and volume need not be 
considered for calculations. 



III. SOME ADDITIONAL PROPERTIES 



We define a function F 



N+l 
3=1 



(14) 



and see from Eq. ([6]) that 

V 2 F = 0, F(Si) = l (i=l,...,N+l). (15) 

Since F — 1 throughout the surface St, we see by 
uniqueness that F = 1 in the volume Vs T from which 
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we find that 



N+l 

i=i 



(16) 



In addition, by summing over j in Eq. (|7b[) and taking 
into account Eq. (fH?)) . we find that 

^C«=0 (i=l,...,iV + l). (17) 

3=1 

The symmetry of the Cij elements leads also to 

AT+l 

£)C«=0 (J = 1 JV + 1). (18) 



Equations (jT7|) and |T5J imply that the sum of the ele- 
ments over any row or column of the matrix is zero. Ap- 
pendix [5] gives some proofs of consistency for these im- 
portant properties. Taking into account the symmetrical 
nature of the Cij matrix with dimensions (N+l) x (N+l) 
and the N+l constraints in Eq. |T5J), we see that for 
a system of N conductors surrounded by another con- 
ductor N + l, the number of independent capacitance 
coefficients is 



(A+l) 2 



N(N + 1) 



Other important properties are that 



C u > 
Cij < 0, 



(i + 3). 



(20a) 
(20b) 



Equation (|20a[) follows straightforwardly from Eq. (JTT 
To demonstrate Eq. (|20b|) , we recall that the solutions 
of Laplace's equation cannot have local minima nor local 
maxima in the volume in which the equation is validi^ 
Consequently, the fj functions must lie in the interval 



< /,■ < 1. 



(21) 



Because fj = on any surface Si for i ^ j, we see that fj 
acquires its minimum value on such surfaces. Therefore 
the function V/j should point outward with respect to 
the conductor i for i ^ j. Hence 



n i • V/j > for i ^ j. 



(22) 



We substitute Eq. ^ into Eq. © and obtain C!y < 
for i ^ j. An additional derivation of the fact that 
Ca > can be obtained by taking into account that fj 
acquires its maximum value on the surface Sj . 
Equation (TTB"]) can be rewritten as 



V a, = -c 



i=i 



N+1J- 



(23) 



From Eq. (f2T)]) we have that Cn+ij < for j = 1, . 
and Cjv+i,jv+i > 0. Hence 



, N 



N 
i=l 



(3 



,N) 



N 



(24a) 



(24b) 



The following properties follow from Eqs. (HJ), P^|) . 
(1201), and 



N 

|c«| > E |Cy 

Cu | > | Cjj | , 
CuCjj > cij 

N 

|CV+l,Ar+l| = 2^ \^i> N + 1 \ 
i=l 

|Cjv+i,iv+i| > |Cj,jv+i|, 



(25a) 

(25b) 
(25c) 

(25d) 

(25e) 



where i,j = l,...,N. 

A particularly interesting case arises when the external 
conductor is at zero potential. In such a case, although 
the elements of the form Cm+i ,j do not necessarily van- 
ish, they do not appear in the contributions to the charge 
on the internal conductors as can be seen from Eq. ([7]) by 
setting (fN+i = 0. For this reason, the capacitance ma- 
trix used to describe N free conductors (that is, not sur- 
rounded by another conductor) has dimensions N x Nr^ 



IV. TWO CONDUCTORS 

We illustrate our method by deriving the basic proper- 
ties of a system of two conductors. These examples will 
show the usefulness of Eq. (|7|) and some of the properties 
derived from our approach. We analyze a single internal 
conductor with an external conductor that is, N = 1. 
The internal conductor is labeled as conductor 1. From 
Eqs. HJ) and ([18]) we have 



C'21 — C\2 — —Cu — —C; 



22- 



(26) 



Therefore, there is only one independent coefficient, say 
Cu (in agreement with Eq. ^ with N = 1). The 
charges on the internal and external conductors can be 
calculated from Eq. ([7]) 



11 = Cn(tpi - tp 2 ) 

12 = -C n (ipi - ¥2) = -Qi- 



(27a) 
(27b) 



Equation (|27b[) is consistent with Eq. (| A4[) and shows 
that the charge induced on the surface of the cavity of the 
conductor 2 is opposite to the charge on the conductor 
1. 
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In Table U we display the results of three well known 
configurations of two conductors. The second column 
shows the fi functions, which can be found by Laplace's 
equation © and used to calculate Cu with Eq. 0. 



System 


fx 


C u 


Spherical shell with radius b and 
concentric solid sphere with radius 
a. 


ab (1 _ 1\ 
b—a V r b J 


47T£q ab 
b—a 


Cylindrical shell with radius b and 
concentric solid cylinder with ra- 
dius a, both with length L. 


]n(r/b) 
hx(a/b) 


2tt£q L 
ln(6/o) 


Two parallel planes with area A at 
x — and x — d (conductor 1). 


x/d 


e A 
£ 0~d 



TABLE I: Cu and fi factors for three systems of two con- 
ductors with a < r < b and < x < d. We neglect edge 
effects for the cylinders and planes. 



V. EXAMPLES 

We use our approach to study a system with embed- 
ded of conductors. In addition, the case of two inter- 
nal conductors is examined, and we show the limit in 
which the configuration of two conductors without ex- 
ternal conductor is obtained. These examples show how 
the properties we have derived can be used to calculate 
the capacitance coefficients. 

Example 1. Consider two concentric spherical shells 
with radii b and c and a solid spherical conductor (con- 
centric with the others) with radius a such that c > b > 
a. The potentials are denoted by tp\, (f2, and (p 3 respec- 
tively. The general solution of Laplace's equation for fi 
can be written as 



r 



B,. 



From Eqs. ([6]) and (|28j) we obtain /i and / 3 




(a < r < b) 
(b<r <c) 



h = 



be 
c-b 



(a < r <b) 
(b<r<c) 



(28) 



(29) 



(30) 



Although fi can be obtained the same way, it is easier 
to extract it from Eq. (TTB]). The result is 



) (b < r < c) 



be tl 
c—b V r 



(31) 



The nine capacitance coefficients can be evaluated ex- 
plicitly from Eq. |(7J), but it is easier to use Eqs. (fTl]) and 



(fTg)) and to take into account that C31 = (V/i(r) 
for r > b). We have 



C13 - 0, 
C22 = Cu 



C: 



:>,2- 



C12 

C33 



-Cu 
-C32. 



From Eq. jJJ the charge on each conductor is 



Qi = Cix{ipx - tp 2 ) 

Q2 = -Qi + (732(^3 - fi) 

Q3 = C 3 2(^2 - Va) = -{Qi 



Q2 



(32a) 
(32b) 



(33a) 
(33b) 
(33c) 



Hence, we only have to calculate Cu and C32*^ The 



result gives 



Cu = 4irs - 



ab 



c 



32 



-4tt£ - 



bc 



(34) 



If (f2 — </?3, we find that Q\ — —Q2 and Q3 = 0. 
It can be shown that Eqs. (|32|) and (f33|) are valid even 
if the conductors arc neither spherical nor concentric, 
because those equations come from Eqs. ([Tail . ([IT]) , and 
(|18J) which are general properties independent of specific 
geometries. 

Example 2. Consider two internal conductors and a 
grounded external conductor. As customary, we begin 
with Qi = Q2 = 0. By transfering charge from one in- 
ternal conductor to the other we keep Q\ = —Q2- From 
Eq. (|7a|) and defining V = tpi — tp2 we find 



(35) 
(36) 



Qi = {Cu + C12) tpi - C12V, 
Qi = -Cnipi - C 12 V, 

where we have used Eq. (jT5J) . Similarly Q2 - 
C22V, and using again Eq. (fTS"]) we find 

Ql + Q2 = C , 33^ 1 -C 3 2K 



-C23P1 



(37) 



Since the system is neutral Q\ + Q2 — and hence 

<Pi = (38) 
L-33 

substituting Eq. (|3^|) into Eq. (|36[) we obtain 
Qi = CV 



C = C ^~ C ^\ (39) 



C33 

Because N = 2 only three of the coefficients in the 
definition of C are independent. From Eqs. ([20|) we 
see that this effective capacitance is non negative. The 
procedure is not valid if C33 = 0, in that case we see 
by using Eqs. (H} and ([20]) that C i3 = C 3i = 0, and 
from Eq. (f3"6) we find C — — C\ 2 — C22, which is also 
non negative. The limit in which there is no external 
conductor is obtained by taking all the dimensions of the 
cavity to infinity while keeping the external conductor 
grounded as discussed in Ref. [ll|. 
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VI. CONCLUSIONS 



We utilize Eq. to write Qi nt as 



We have used an approach based on Laplace's equa- 
tion to demonstrate that the capacitance matrix depends 
only on purely geometrical factors. The explicit use 
of Laplace's equation permits us to demonstrate many 
properties of the capacitance coefficients. The geomet- 
rical relations and properties shown here permits us 
to simplify many calculations of the capacitance coeffi- 
cients. We emphasize that Laplace's equations necessary 
for finding the capacitance coefficients are purely geo- 
metrical as can be seen from Eqs. ^ and 0. Laplace's 
equation is usually easier than Green function formalism 
for either analytical or numerical calculations. Appendix 
[X] shows some proofs of consistency to enhance the phys- 
ical insight and the reliability of our method. 



APPENDIX A: PROOFS OF CONSISTENCY 

A proof of consistency for the identity (fTg)) , is achieved 
by using Eq. to calculate the total charge on the N 
internal conductors^ 



N+l 



Qint — ^ Qi — 



N 



We use Eq. (118]) to find 

Qint = _ 



3=1 L i=l 



N+l 



(Al) 



(A2) 



Note that Eq. (|A2[) requires many fewer elements of the 
Cij matrix than Eq. (IA1[) . This difference becomes more 
significant as N increases. If we again use Eq. 0, we can 
find the charge on the cavity of the external conductor 



N+l 

}n+1 = CN+l,j<Pj, 
3 = 1 



and therefore 



<3jv+i 



(A3) 



(A4) 



a property that can also be obtained from Gauss's lawi^ 
Another proof of consistency for Eq. (fTS)) is found by 
employing Eqs. and (|A2|) to calculate Qint (taking 
into account that Eq. (|A2[) comes directly from Eq. (fT8|) ) 




(A5a) 

n N+1 dS. (A5b) 



Qint=e j> \7<j)-n N+ i dS = e j> E-(-n N+1 ) dS. 

J Sn+i J Sn+i 

(A6) 

This relation is clearly correct because njv+i points in- 
ward with respect to the volume Vs T ■ 

A proof of consistency for Eq. p0[) that shows the sym- 
metry of can be obtained by calculating the electro- 
static internal energy, which in terms of the electric field 
is 



U 



so 



f E 2 dV= 6 -^ f X7cj)-V<j)dV (A7a) 

Jv St 2 Jv St 



IV Sl 
N+l 



- V tpupj e / V/i • V/j 
1 : ' Jv.^ 



dV 



(A7b) 



where we have used Eq. ([5]). From Eq. (JTUJ) we find 



N+l 



N+l 



U = \ C ij<P3<Pi = \ Q 



(A8) 



h3 



consistent with standard results- 1,2 



APPENDIX B: SUGGESTED PROBLEMS 

To enhance the understanding of this approach and 
its advantages, we give some general suggestions for the 
reader. 

1. Implement a numerical method to solve the 
Laplace's equation ^ for the functions associ- 
ated with a nontrivial geometry (for example, two 
non-concentric ellipsoids). Use Eqs. (fTI))) and (fJT|) 
to either simplify your calculations or to check the 
consistency of your results. Then use Eq. to ob- 
tain the Cij factors numerically. Use Eq. (jlip and 
Eqs. (fT7l) -([25 |) either to simplify your calculations 
or to check the consistency of your results. 

2. We have emphasized that to calculate the total 
charge on the internal conductors Eq. (|A2j) requires 



many fewer Cy elements than Eq. (|A1 [) . How 
many fewer elements are required for an arbitrary 
value of N? 

3. For a successive embedding of concentric spherical 
shells, calculate the capacitance coefficients for an 
arbitrary number of spheres. 

4. Show that for the successive embedding of three 
conductors with arbitrary shapes, Eqs. (|32p and 
([33]) still hold. Generalize your results for an arbi- 
trary number of conductors. 
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Qn+i is not necessarily the total charge on the external 
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obtain an integral of the volume directly from Eq. (0. 
However, fj is not defined in the region inside the con- 
ductors. The gradient of fj in Eq. is evaluated in an 
external neighborhood of the conductor surface. 

11 By uniqueness, the solution for this problem is equiva- 
lent to the solution for a system consisting of the same N 
conductors contained in the cavity of a surrounding con- 
ductor, such that all the dimensions of the cavity tend to 
infinity, and the potential of the external conductor is set 
to zero. 
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dom (although zero), we have a total of three degrees of 
freedom, in agreement with Eq. (|19|l for N — 2. 

13 For a derivation of some of these results based on the en- 
ergy of the electrostatic field see L. D. Landau, E. M. Lif- 
shitz, and L. P. Pitaevskii, Electrodynamics of Continuous 
Media (Elsevier Butterworth-Heinemann, 1984), 2nd ed., 
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